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Abstract 

The problem of weakly correlated electrons on a square lattice is studied theoretically. A simple 
renormalization group scheme for the angle-resolved weight Z(ff) of the quasiparticles at the Fermi 
surface is presented and applied to the Hubbard model. Upon reduction of the cutoff the Fermi 
surface is progressively destroyed from the van Hove points toward the zone diagonals. Due to 
the renormalized Z(9), divergences of both antiferromagnetic and superconducting correlation 
functions are suppressed at the critical scale, where the interactions diverge. 
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Understanding of the one-particle spectrum of strongly correlated systems near the 
metal-insulator transition is an extremely difficult task, particularly if one wants to con- 
struct a microscopic theory. One standard example is the pseudogap regime of the HTC 
superconductors. ARPES measurements |IJ showed that the Fermi surface is destroyed by 
correlations. This happens first near the van Hove points where the one-particle spectrum 
develops the characteristic 2-peak structure. Remaining parts of the Fermi surface, often 
called Fermi patches, get progressively narrower around Brillouin zone diagonals as the 
temperature decreases. Regions around van Hove points contain non-Fermi liquid with a 
pseudogap and other signatures of strong correlations such as flat bands. In other words, 
the pseudogap has the form similar to the absolute value of the ^.^-superconducting (SC) 
order-parameter. This vision of the pseudogap is in agreement with STM results as well, 
regardless of details on how these results are interpreted. The above experiments still do 
not reveal much on the origins of the pseudogap, namely whether it is simply the signature 
of a liquid of pre-formed pairs or something much richer in fluctuations. In fact it is known 
that the antiferromagnetic (AF) fluctuations are also strong in the pseudogap regime. [f§] 

If we assume that AF and SC fluctuations are somehow together the major reason for 
the strong renormalization of the one-particle selfenergy, and for the consequent partial 
destruction of the Fermi surface, then a many-body analysis of the one-particle propagator 
can be done in a controlled way. In fact the weak coupling theory easily reproduces AF and 
SC fluctuations from particle-hole (p-h) and particle-particle (p-p) loop-logarithms. Even 
if the coupling in realistic HTC systems is of the order of Fermi energy (i.e. intermediate-to- 
strong), already the weak-coupling theory contains the observed two-particle correlations. 
In this paper I answer the question of how the angle-resolved quasiparticle weight Z{9) is 
renormalized by strong and coupled AF and SC fluctuations, and of the main consequences of 
the renormalized Z(8) to the characteristic angle-resolved two-particle correlation functions. 
The renormalization of the quasiparticle weight in 2D was recently studied by Kishine and 
Yonemitsu. || To calculate the renormalization of Z resolved in the position on a flat Fermi 
surface they used a two-loop selfenergy expansion with the two-loop-renormalized vertex. 
The results show clearly that the flatness of the Fermi surface induces the suppression of the 
quasiparticle residue and that this effect is anisotropic. In the present work I consider the 
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whole square Fermi surface of the Hubbard model. For this purpose I employ the N-patch 
renormalization group theory. 

After recent theoretical studies from several groups it emerged that the N-patch model 
describes in a systematic and controlled way weakly correlated electrons near half-filling, 
and explains the major aspects of HTC-s.|| |7|, |[ |], [10], |ll| Until now the RG analysis of the 
N-patch model has been done only on the level of the two-particle scattering amplitudes 
or, in field-theoretical jargon, of the four-point vertex U(Ki, K2, K3). The analysis of the 
renormalization group flow of U as a function of three patch indices ii,*2)*3 gave several 
important results. Typically the amplitudes ?7(zi, 22, 23) diverge at some interaction- and 
doping-dependent critical energy scale A c . For the case of the Hubbard model, we distinguish 
two main renormalization regimes ||, the parquet regime and the BCS regime. 

In the parquet regime (|/i| < A) both particle-particle and particle-hole propagators 
have strong contributions to the beta-function due to the van Hove singularities and the 
Fermi surface nesting. In this regime, provided A — > A c and neglecting the selfenrgy cor- 
rections, both SC and AF tendencies are strong and build divergent correlation functions 
X and X ■ The dominant component of the antiferromagnetic susceptibility is of the 
s-type and the dominant component of the superconducting one is of the rf a .2_ y 2-type. Both 
static compressibility Xc and homogeneous magnetic susceptibility Xs go to zero as the cut- 
off approaches to its critical value A c .|| [| |TU| Consequently A c is the energy scale of the 
crossover between the strange metal and the strongly correlated regime with gap or pseu- 
dogap. A question arising from already existing results on the Hubbard model || [8] and 



on its extensions || [II], |TT| is the following: if one wants to interpret the critical scale A c 
(or temperature) as the energy (temperature) T* for the onset of the pseudogap, why then 
are all signatures of the pseudogap not seen? This means in particular that x AF an d X SC 
should be finite and not diverging. Emergent is the necessity to calculate the correlation 
functions with the corrections due to the one-particle selfenergy. 

At stronger doping (\fi\ > A) nesting properties get weaker so that eventually only re- 
maining renormalization channel is superconducting (p-p). This is the BCS regime. There 
the superconductivity is simply BCS-like with the coupling constants and the angular pro- 
file of the order parameter determined at higher scales by a parquet-like flow, where A was 
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larger than the chemical potential. 



Equivalent to the RG approach is the fast parquet theory |12], [13| , where 6 variable 
(continuum version of the patch index) is called the fast variable in addition to the cut- 
off logarithm called the slow variable. In two dimensions the parquet integro-differential 
equations always have mobile pole solutions, i.e. the AF and SC fluctuations decouple one 
from the other. This seems to be in disagreement with the RG results, where we detected 
only immobile poles, the type of solution in which all scattering amplitudes develop the 
pole at the same scale A c . The question of the consilience between the two theories is still 



controversial. However, the results of De Abreu and Dougot |TJ| indicate that the mobile 
pole solution is dominant only in the very vicinity of A c and that the fixed pole solution 
is an intermediate solution, valid over several decades of energy scale. The width of the 
"very vicinity" characterized by the mobile poles depends on the coupling constant so that 
for reasonable and not too weak Uq the final regime is so close to A c , that the couplings 
are already too strong and out of reach of a weak coupling theory. Consequently the real 
physical interpretation can be given only to the immobile pole regime. We will concentrate 
on this "intermediate" regime in which all fluctuations are coupled and at least behave as if 
having an immobile pole. 

We suppose that the electronic Green function has the form 

«w - f (k) 

Zi(9) is the angle-resolved scale dependent quasiparticle weight and £(k) is the tight-binding 
dispersion. The formalism keeps the notation introduced in the reference ||. The form ([!]) 
contains two main approximations. The first one is to keep trace only of the renormalization 
of the coherent part of the propagator. The second approximation is to assume that the 
spectrum £(k) remains non-renormalized. This assumption implies that we ignore the flow 
of the Fermi surface (FS) and of the Fermi velocity. Because of the particle-hole symmetry 
the flow of the FS is zero at half-filling, where we can expect that our form of the Green 
function is closer to reality than in the imperfectly nested (non-half filled) case. 

The flow equation for Ziff) is derived from the general and exact one-loop RG equation 
for the complete selfenergy S(K) given in M. Let us suppose that we are at some scale I, 



and that we know the propagator ([I]) . We integrate dl further and look what is the effective 
two-point vertex T 2 in the effective action S(l + dl); it is 

T 2 (l + dl) = Z- 1 (6)(ilu - £(k)) + dE t (K) . (2) 
To find Zi+di we expand dHi{K) in first order in iuj to obtain 

Z l+<B = Z l (l-Z l d iu dZ). (3) 
The differential equation for Z follows immediately 

diZ x {Q) = -Zf(0) [$£(#)] \^ =0 . (4) 
Only the terms of diH which are linear in energy contribute. These are just the terms which 



are marginal upon zeroth order scaling in Shankar's sense. [|T|] We will look for these terms. 
The equation for <9;£ can be written as 

d&{0,e,u>) = ^i/^E/ M0',A)d9'G l (9',co',uA)V l (K,K' u ) , (5) 

where T>\ = 2Fi — Ff, Fi and F] are energy-momenta dependent forward and back- 
ward scattering processes at the scale I, related to the effective interaction in a way that 
Fi(K lt K 2 ) = U{K U K 2 , K x ) and F^i, K 2 ) = U(K 1 , K 2 , K 2 ). J v {6, A) is the angle-resolved 
density of states at energy £ = vk. There is another, approximate but physically justified 
way to decompose T>. In fact, we will suppose that Di can be written as a sum of p-p and 
p-h terms: 

V t (K, K') = Vf[K + K') + Vf{K - K') . (6) 

The p-p part of the propagator T> depends only on the total energy-momentum Q pp = 
(cjpp, q pp ) = K + K' while the p-h part depends only on the energy-momentum transfer 
Q P h = (u p h, q p h) = K — K' . As usual we skip the marginal part of the dependence on q pp 
and (\ p h. For that purpose we note that both momenta can be written in the form 

q = q (0) (M0+q (1) (M\£,O, (7) 

where q stands either for q pp or for q^, c^°\9, 9') is the value of q when both momenta k 
and k' are at the Fermi surface, while q^ 1 -* is the correction due to non-zero energies £ and 



Using standard scaling arguments we can skip in the limit A/ep — > because the 
cutoff is imposed to momenta. Similar argument cannot be used for u pp and u p h because 
the integral in (|5|) runs over all frequencies independently of the actual cutoff A(l). We are 
therefore left with the locally dispersionless phononic propagators 

Vf p (6, 9', Q pp ) « Vf(9, 9', iu + too') = 2F PP (9, 9', iu + iu') - F pp (9, 9', iu + iu') (8) 

and 

V ph (9, 9\ Q ph ) « V ph {9, 9', iu - iu') = 2F ph (6, 6', too - too') - F[ h (9 7 6', iu - iu') . (9) 

In these expressions we made the same pp-ph decomposition of the forward and backward 
amplitudes as we did with V in eq.||. The following step is to re-constitute the iu dependence 
from the cutoff dependence. This can be done with logarithmic precision simply replacing A 
with iu. The derivatives over frequency in eq.(|J) (acting only on F-parts) are then readily 
calculated 

d i uF( p * h {e,e\iu±iu')\ iw= v = ±^Ff^(M') , (io) 

and equivalently for backward amplitudes Ff p,ph . The frequency-independent quantities 
diF[ p (8, 6') and diFf h (6, 9') are the p-p and p-h parts of the /3-function of the N-patch 
model, with appropriate configurations of the external momenta: 

d l F pp (9,9')=P pp {U,U}(9,9',9) 
8^(9,9') =[3 pp {U,U}(9,9',9>) 
8^(9,9') = -[X(3 ph {XU,XU}}(6,6',0) 

dtF ph (9,e') = [2(3 ph {U,U} -f3 ph {U,XU} -f3 ph {XU,U}}(9,9',9') , 
where all /3-functions are given in ref.||, but for the moment with dressed Green functions 
(|IJ). Notice that the forward scattering has finite p-h contributions only from the ZS' channel 
(1=3) while only the ZS channel (1=4) contributes to the backward scattering. This means 
that we forget about the contributions at zero momentum transfer. They are somewhat 
tricky, but don't have any logarithmic part so that we can forget them. 

We can also get rid of the Z factors in beta-functions of the eq.flTTD by rescaling the 
fermions at every step of the RG in a way that 

^{9)Z~ 1 {9)^{9) -> *(0)tf (0) (12) 
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and re-defining the effective interaction 

E/,(l, 2, 3) -> [Z i (l)^(2)^(3)^(4)]- 1 / 2 ^(1, 2, 3) . (13) 

After transformations (0) and (|P2| ) the calculations of the /3-functions to the one-loop order 
are identical to the case with Z — 1. Performing uj' integral to logarithmic precision, the 
flow equation for Z{9) becomes 

d l \ogZ l {e) = J d8' J_(8',A)rii(9,9') =T}i(9) , (14) 

where 

Vl (9, 9') = di{2[Ff - Ff h ] - Ff + F[ h }(9, 9') 

with 9;F-terms given by eqs. flTl|) and calculated with the bare Green functions, as in ref. ||. 
Generalization of eq.[l4] to finite temperatures can be done simply by replacing J71 with 
E„(-^^njr(vA). 

Taking ID "limit" of the eq.(|T4]) is simple and instructive: instead of N patches we now 
have 2 patches: 9 = R (right) and 9 = L (left). Independent scattering amplitudes at 
non-rational filling are F(RL) — g 2 — U(RLR) and F(RL) — g± — u(RLL). We skip 
#4 = u(RRR) from considerations because it has no logarithmic renormalization. The 9- 
integrals reduce to summation over two points so that it is easy to reproduce the well-known 
result T] = — I 2 (g\ — gig 2 + gT) |l^| . It is the Luttinger liquid exponent. 

In two dimensions the angle resolved 77^ (6*) can also be associated with some non-Landau 
(non-Fermi liquid) behavior. We see that r\ becomes finite if the forward and backward 
amplitudes have some logarithmic flow over a wide range of (9-space (A9 ~ 1). Another 
possibility for having finite rj is when the Fermi surface is close to the van Hove singularities. 
"Close" means that the distance between the Fermi level and the van Hove singularity is 
comparable or inferior to the scale A c at which interactions start to flow strongly. 

We will now calculate the renormalization of Zi{9) in the 2D Hubbard model at half- 
filling, from the knowledge of the scale dependence of the patch-dependent interaction 
Z7j(l,2, 3). In the RG equations of the previous section the discretization of 9 is done in 
a way described in ref.||. In the present case the Fermi surface is square so that there are 
two mechanisms for the suppression of quasiparticle residues. Namely, both above mentioned 
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conditions are fulfilled: (i) forward and backward amplitudes have logarithmic flows for any 
configuration (9, 9') if the two angles are at opposite sides of the Fermi surface, so that the 
available phase space is indeed large; (ii) van Hove singularities are at the Fermi surface and 
are nested. In fact, one can alternatively imagine a Fermi surface with non- nested van Hove 
singularities and nested parts elsewhere. Such a model would be even closer to the realistic 
situation in some HTC compounds. For the sake of rigor, we will however remain limited 
to the Hubbard model. 




FIG. 1: The evolution of the angle-resolved quasiparticle weight on the Fermi surface. The lines 
are for I = log (At/ A) =3.; 4.; 4.4; 4.95; 5.11; 5.17; 5.20. The critical scale is l c « 5.204 

The result is shown in fig.|l|. The figure shows Z(9i) with 0<#j<7r/2on9 equidistant 
points. Different lines correspond to different values of the scaling parameter as it approaches 
its critical point, i.e. when couplings diverge. Settings are the same as in ref.||: the Fermi 
surface was discretized into 32 patches and the initial interaction is Uo/(4t) = 0.333. One 
sees that the Fermi surface is first destroyed at the van Hove points and than the regions of 
the FS destroyed by correlations grow larger and larger. This kind of flow is compatible with 
the interpretation that A c is not the critical temperature for some symmetry braking, but 
merely the scale at which coherent quasiparticle cease to exist at the Fermi surface giving 
place to a gapped or pseudogapped liquid. The magnitude of the (pseudo-)gap is largest in 
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van Hove points and smallest on the diagonals of the Brillouin zone. The whole gap function 
A(0) then scales as A c x f(9), where f(9) is a function with the symmetry of the absolute 
value of the d x 2_ y 2 -harmonic. This is the angle-resolved (pseudo-)gap, responsible for the 
correlation-induced angle-resolved localization. In other words the electrons near the van 
Hove points get much less mobile than those near diagonals. The similar scenario has been 
proposed by the Zurich group [TO] even without concrete calculations of the quasiparticle 



weight. 

The question of antiferromagnetism and superconductivity remains to be clarified. Let's 
discuss this problem taking into account the scale dependent Z(9) in the flow equations for 
the susceptibilities x AF (9,9') and x SC (9,9'). Following the procedure given in ref.|| and 
dressing the electronic propagators with Z-factors we get 

*f(*l, fc) = Zl{ QJ Zl{ Q 2) j dd 0)afP)3f(0. AO • ( 15 ) 

This equation has the same structure as the one in reffl, with two modifications. First, we 
skip the retardation effects, replacing U simply by I, because we are at half-filling. Second, 
the quantity zf(9i, 9) that has the role of a triangular vertex is somewhat modified. Its flow 
writes: 

- V (9i) - 77(0 2 )] 9 2 ) = -j d9 zf(9 1 , 9)D*(9)V l s (9, 9 2 ) . (16) 
The meaning of zf(9i, 9 2 ) is that 

zf(9 1 ,9 2 ) = Z l (9 1 )zf(9 1 ,9 2 )Z l (9 2 ) 
so that the initial conditions for z\ and for z\ are the same. After discretization we integrate 



numerically equations ([TBI ) an d (16). Fig|J shows the flow of the dominant eigenvalues of 
susceptibilities 5 = AF and 5 = SC near the divergence of scattering amplitudes. The thin 
line represents both (degenerated) susceptibilities for U = 0. Including only the one-loop 
vertex renormalization we get the strong enhancement and, as far as my numerics can say, 
even divergences of both AF and SC susceptibilities. If we now include also the one-particle- 
weight renormalization, both susceptibilities are radically reduced and lose their divergent 
behavior. On the other hand, the flow of the compressibility Xc and magnetic susceptibility 



Xo- is not affected by Z-renormalization because of the Ward identities. |17| All above results 
support the statement that what happens at energy scale A c is a flow towards a state with 
spin- and charge-gap or pseudogap, insulating and without AF or SC ordering. 




FIG. 2: Scale dependence of the dominant components of both antiferromagnetic and superconduct- 
ing susceptibilities at half-filling, A: due to renormalized vertex only and B: due to renormalized 
vertex and quasiparticle weight. The thin line represents the bare susceptibility. 

We finish this discussion with a few words about the effects of doping. Two regimes, 
parquet and BCS, exist also in the scaling properties of Z\{&). In the whole parquet regime 
we expect the behavior governed by the proximity to the half-filling situation, so that the 
present results can be applied at all energies larger than the chemical potential \\x\. Upon 
doping the nesting becomes more and more imperfect, the p-h logarithm loses its divergence, 
the critical scale is more and more suppressed and the p-h and p-p channels get progressively 
less coupled. Eventually at strong enough doping and low enough A one is in the BCS regime 
where the effective physics is described by the 2D BCS theory, with renormalized and 6- 
dependent interaction and quasiparticle weight. The anomalous dimension rj(9) goes to zero 
because (i) the the range Ad, over which the forward and the backward amplitude have 
strong flow due to BCS diagram, scales with the cutoff, and (ii) the van Hove singularities 
are outside of the cutoff. The critical temperature for the onset of the superconductivity 
is not affected by the renormalized Z, but the magnitude and angular dependence of the 
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superconducting order parameter are dressed by Z(9). The symmetry of the gap remains 

d x 2 —y 2 • 

To summarize, I proposed a simple renormalization group theory for the angle-dependent 
destruction of the Fermi surface in the Hubbard model. The results offer a theoretical com- 
prehension of the angle-dependent Fermi surface truncation in the cuprate superconductors 
in terms of the scattering processes of the electrons on the low-energy collective excitations 
of both particle-particle and particle-hole types. The theory, based on the N-patch model, 
is in its essence a controlled weak-coupling procedure that keeps trace of the dependence 
of the effective interaction and one-particle spectral weight on the position of the particles 
at the Fermi surface. As one approaches the critical scale A c the quasiparticle weight goes 
to zero first near the van Hove points, and the effect progresses toward Brillouin zone diag- 
onals as one lowers the temperature. Dressing the flow equations for AF and SC response 
functions with the one-particle weight factors results in dramatical reduction of correlations 
of both types. The strongly correlated state just below A c is gapped or pseudogapped and 
without any long-range-order. Critical scale A c is interpreted as the pseudogap temperature 
T* found in cuprate superconductors. 
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